Simulation of point defect diffusion in structures with local elastic stresses 
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The stress-mediated diffusion of nonequilibrium point defects from the surface to the bulk of 
the semiconductor is investigated by computer simulation. It is supposed that point defects are 
generated in the surface region by ion implantation and during diffusion pass over the local region 
of elastic stresses because the average defect migration length is greater than the thickness and 
depth of the strained layer. Within the strained layer point defect segregation or heavily defect 
depletion occur if defect drift under stresses is directed respectively in or out of the layer. On the 
other hand, the calculations show that, in contrast to the case of local defect sink, the local region 
of elastic stresses practically does not change the distribution of defects beyond this region if there 
is no generation/absorption of point defects within the strained layer. 
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I. INTRODUCTION 



II. MODEL OF STRESS-MEDIATED 
DIFFUSION OF POINT DEFECTS 



In recent years decreasing dimensions of the inte- 
grated components especially with various multi- 
layered structures^ such as Si/SiGe ; 4 i 5 i 6 and the use 
of different nano-particles embedded in the crystalline 
matrisZi^, have been the main trends of advanced de- 
vice technologies. Much attention is given to studying 
of the defect kinetics and stress evolution during semi- 
conductor processing because defect and stress engineer- 
ing can significantly improve the device performance ^2*2 
For example, many efforts are directed to a study of 
the stress-mediated diffusion of impurity atoms and 
point defect Sfifi^i*iSiiii4 including the influence of elas- 
tic stresses on the drift of point defects near the sur- 
face or at the interfacesji^iSiii and near different inho- 
mogeneities of semiconductor crystalsiifi^. Also, many 
studies are associated with changes in the defect sub- 
system of ion-implanted layers as these changes are re- 
sponsible for the transient enhanced diffusion of dopant 
atoms (iiSiiiLSfli For example, to explain the experimental 
data, it is assumed^ that elastic stresses arising in the re- 
gion of the SiGe layer buried in silicon and compressively 
strained after growth cause the drift of vacancies to this 
layer. Thus, accumulation of the vacancies within the 
layer and their transformation into nanovoids occur. In 
RefiiS the transient enhanced diffusion of dopant atoms 
in hyperfine boron-doped layers created by the molecular- 
beam epitaxy is investigated experimentally. It is as- 
sumed that the formation of the clusters of boron atoms 
with silicon self-interstitials occurs in these doped regions 
during thermal treatment. This process and a number of 
similar thermal treatments are characterized by the local 
absorption of self-interstitials. Thus, investigations con- 
cerning the influence of the local strained regions and lo- 
cal sinks on the diffusion of nonequilibrium point defects 
are of great importance for the next-generation fabrica- 
tion processes. 



To calculate distributions of the nonequilibrium point 
defect concentration in the field of elastic stresses, the 
stationary diffusion equation established in Refi^ can be 
used. In the form convenient for numerical solution this 
equation can be written as^i 
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Here C x and Cf are the concentration of point defects 
(vacancies or self-interstitials) in the neutral charge state 
and equilibrium value of this concentration, respectively; 
d(x) and di are the effective diffusion coefficient of point 
defects and intrinsic diffusivity of these defects, respec- 
tively; C is the concentration of substitutional^ dissolved 
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dopant atoms forming a doped region; Cb is the total 
concentration of dopant atoms responsible for opposite- 
type conductivity; ni is the intrinsic carrier concentra- 
tion; k c (x) is the normalized to fej concentration depen- 
dence of the effective absorption coefficient for point de- 
fects; k Sp (x) is the spatial distribution of the absorption 
coefficient; ki and T{ = k^ 1 are the values of absorption 
coefficient and average lifetime for point defects in intrin- 
sic semiconductor, respectively; U = \fd^ji is the average 
migration length of point defects in intrinsic semicon- 
ductor; v x is the x-coordinate projection of the effective 
drift velocity of point defects due to elastic stresses; U d is 
the potential energy of these defects in the field of elas- 
tic stresses; g R is the rate of the noncquilibrium point 
defect generation per unit volume due to the external ra- 
diation or dissolution of extended defects; gi is the rate 
of thermally equilibrium generation of point defect in in- 
trinsic semiconductor. The function a d — h r /h\ takes 
into account that the real constants h r for the transition 
between the charge states of point defects deviate from 
their equilibrium values h\ due to heavy doping effects 
and elastic stresses. 

Eq. ||TJ differs from the defect diffusion equation used 
in Reff2S because the drift of mobile species in the field of 
elastic stresses is included. On the other hand, from Eq. 
l|Tjl the distributions of nonequilibrium point defects can 
be derived, instead of computing the equilibrium defect 
distributions by means of expressions used in Ref>i2iii 
Eq. (yi is also convenient for numerical solutions due 
to the following characteristic features: (i) This equa- 
tion describes the diffusion-reaction-drift of nonequilib- 
rium point defects with different charge states as a whole, 
although only the normalized concentration of the neu- 
tral defects C must be obtained to solve the equation, 
(ii) The obtained equation takes into account the drift 
of all charged species due to the built-in electric field. 
At the same time, there is no explicit term describing 
this drift, (iii) Despite the fact that the effective coeffi- 
cients of Eq.Q represent nonlinear functions of Xi these 
functions are smooth and monotone 

It is also assumed that the mobility of point defects 
is significantly greater than the impurity atom mobility, 
and there is no change in the processing conditions or 
the changes are sufficiently slow. In this case the time 
derivative of the defect concentration is close to zero and 
distributions of point defects are quasi-stationary with 
respect to the distributions of dopant atoms, clusters, 
extended defects, and also with respect to changes in the 
processing conditions. 



III. NUMERICAL SOLUTION 

The finite-difference method^ 2 , is used to find a nu- 
merical solution for Eq. JIJ in the one-dimensional (ID) 
domain [0, xg]. Following RefiSS, the first term in the 
left-hand side of Eq. is approximated by a symmetric 
difference operator of the second order accuracy on the 



space variable x. At the same time, the second term is 
approximated with the first order accuracy by the asym- 
metric forward/backward-difference operator depending 
on the drift direction. On the other hand, if the defect 
flux due to elastic stresses is comparable to or below the 
defect flux due to the concentration gradient, the sec- 
ond term is also approximated by a symmetric difference 
operator of the second order accuracy. Comparison with 
exact analytical solutions for the particular cases of point 
defect diffusion and calculations on the meshes with dif- 
ferent step sizes were carried out to verify the approxi- 
mate numerical solution. 
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FIG. 1: Comparison of the numerical (solid line) and analyt- 
ical (dots) solutions for the equation of point defect diffusion. 



For example, in Fig. ^ th e numerical solution for Eq. 
||TJ in case of the constant diffusivity and constant coef- 
ficient of defect absorption is presented. For comparison 
with the analytical solution^ Gaussian distribution 
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is used to describe the generation rate profile. Here g m 
is a maximum generation rate of point defect during ion 
implantation normalized to the thermal generation rate 
g%'i Rpd and Ai?p d are the position of a maximum of de- 
fect generation distribution and dispersion of this distri- 
bution, respectively. It is supposed that the defect gen- 
eration occurs due to implantation of hydrogen ions at 
energy 20 keV (R pd = 0.198 fjaa and AR pd = 0.0802 ^m 
are taken from Refi^i). Numerical computations are car- 
ried out on the simulation domain [0, Xb], where Xb and 
mesh point number is are equal to 4.0 /itm and 81, re- 
spectively. To obtain a numerical solution, the Dirichlet 
boundary conditions are imposed on C 



C(0) = 1 , 



C{x B ) = 1 
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The function Cth{x) — 1 is added to the analytical 
solution 2 - to satisfy boundary conditions and take 
into account the thermal generation of point defects. As 
can be seen from Fig. ^ the distribution of nonequilib- 
rium point defects obtained by numerical computations 
agrees with the analytical solution proposed in Ref^. 



IV. SIMULATION OF POINT DEFECT 
DIFFUSION 

In Fig. [21 the calculated distribution of nonequilibrium 
point defects in the structure with the local stress field 
providing the drift of defects into the strained region is 
demonstrated. For comparison, the solution of diffusion 
equation Q in case of zero stresses is also shown by the 
dotted line. 
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FIG. 2: Calculated distribution of the concentration of neu- 
tral point defects passing over the region of local stresses (drift 
velocity directed in of the strained layer). A solution for the 
same diffusion equation in case of zero stresses is given for 
comparison. 
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FIG. 3: Spatial distribution of the normalized drift velocity 
of point defects used in calculations shown in Fig. [5] 



To meet the experiments with the transient enhanced 
diffusion, in all the calculations presented here it is as- 
sumed that the generation of nonequilibrium point de- 
fects occurs due to silicon ion implantation at an energy 
of 60 keV (Rp = 0.081 /jm and AR p = 0.033 ^m are also 
from Ref^), i.e. the generation occurs near the surface 
of the semiconductor. Here R p and AR p are the aver- 
age projective range of silicon ions and straggling of the 
projective range, respectively. To describe the defect dif- 
fusion, the average migration length of point defects U is 
chosen equal to 0.7 /*m to be greater than the thickness of 
the strained layer. The reflecting condition for point de- 
fects on the surface and Dirichlet condition in the bulk of 
semiconductor are used in all cases under consideration. 
The distribution of the normalized drift velocity of point 
defects in the field of elastic stresses for the case of defect 
segregation within the local strained region (see Fig. 0) 
is given in Fig. As can be seen from Fig. [31 the lo- 
cal field of elastic stresses takes place between the defect 
generation region and the bulk of the semiconductor. 

Fig. 0] presents the calculated distribution of nonequi- 
librium point defects in the structure with the local re- 
gion of stresses preventing the defect diffusion into the 
strained layer. The distribution of the normalized drift 
velocity of point defects in the field of elastic stresses used 
for this calculation is shown in Fig. [SJ As can be seen 
from Fig. [5J it is supposed that the local stress field also 
occupies a position between the defect generation region 
and the bulk. 
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FIG. 4: Calculated distribution of the concentration of neu- 
tral point defects passing over the region of local stresses (drift 
velocity directed out of the strained layer). A solution for the 
same diffusion equation in case of zero stresses is given for 
comparison. 

As seen from Figs. and the presence of the local 
stresses results either in enrichment or depletion of the 
stress region by point defects. On the other hand, the 
distributions of nonequilibrium defects beyond the region 
of the local stresses are practically unchanged, regardless 
of the stress barrier that should be overcome by the point 
defects migrating in the bulk of the semiconductor. 

A qualitatively different situation takes place in the 
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FIG. 5: Spatial distribution of the normalized drift velocity 
of point defects used in calculations shown in Fig. 



case of strong local sinks of point defects, for example, 
in the case of the silicon structure with an epitaxially 
grown Sii_ x C x layer iSS This can be seen in Fig. show- 
ing the point defect concentration profile calculated for 
the local sink position at the same place as the position of 
local stresses in the previous calculations. Fig. 0demon- 
strates the spatial distribution for the effective absorption 
coefficient of point defects used in the last calculation. 



interesting to note that the concentrations of nonequilib- 
rium defects beyond the region of the local stresses can 
be significantly changed if the generation/absorption of 
point defects takes place in this region. To illustrate, in 
our previous worki^ it is shown that the generation of sil- 
icon interstitial atoms in the region of the local stresses 
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FIG. 7: Spatial distribution of the effective coefficient of point 
defect absorption used in calculations presented on Fig. [(J 



providing the drift of point defects from the surface re- 
sults in supersaturation of self-interstitials in the bulk of 
the semiconductor. 
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FIG. 6: Calculated distribution of the concentration of neu- 
tral point defects passing over the region of local sink. For 
comparison a solution of the diffusion equation for the ther- 
mally equilibrium uniform absorption of point defects is also 
presented. 

A significant decrease in defect concentrations in the 
regions before and after the sink is easily predictable and 
agrees with the experimental data obtained in Refill It is 



V. CONCLUSIONS 

The influence of local elastic stresses on the formation 
of point defect distributions during diffusion of nonequi- 
librium defects from the surface into the bulk of the semi- 
conductor has been investigated by computer simulation. 
Such local regions of stresses can be formed in the mul- 
tilayered heterostructures or semiconductor structures 
with heavily doped layers. Thus, point defects diffusing 
into the bulk must pass over the region of stresses if the 
average defect migration length is greater than the thick- 
ness of the strained layer. Within the strained layer, the 
point defect segregation or heavy defect depletion occur 
if the defect drift under stresses is directed respectively 
in or out of the layer. The numerical calculations show 
that, in contrast to the case of local defect sink, the local 
region of elastic stresses practically does not change the 
distribution of defects beyond this region provided that 
there is no generation/absorption of point defects within 
the strained layers. 
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